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Abstract 

The compressible Xavier-Stokes equations have been used to compute leading 
edge receptivity of boundarv layers over parabolic cylinders. Natural receptivity 
at the leading edge was simulated and Tollmien-Schlichting waves were observed to 
develop in response to an acoustic disturbance, applied through the farfield boundary 
conditions. 

To facilitate comparison with previous work, all computations were carried out at 
ACc, =0.3. The spatial and temporal behavior of the flowfields are calculated through 
the use of finite volume algorithms and Runge-Kutta integration. The results are 
dominated by strong decay of the Tollmien-Schlichting wave due to the presence 
of the mean flow favorable pressure gradient. The effects ot numerical dissipation, 
forcing frequency, and nose radius are studied. The Strouhal number is shown to 
have the greatest effect on the unsteady results. 
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Introduction 


The aerospace design process is complicated by the wide range of factors which 
affect the flight vehicle. Each of these factors, which include propulsion, structures, 
and aerodynamics, must be considered at each step in the design process. Recent 
research on the National Aerospace Plane (NASP) project has shown a very high 
degree of interdependence among these varied design factors. The aerodynamic shape 
of the vehicle forebody, for example, will act as an inlet compression structure for 
the propulsion system. Therefore, the importance of a thorough understanding of 
the flowfield, even though traditionally essential, has now reached an unprecedented 
priority. 

One of the key aspects of the flow is its transition from a stable laminar state to 
one of turbulence. Important parameters such as drag, skin friction, and heat transfer 
differ for laminar and turbulent flows. Small changes in the transition location can 
have significant effects on the overall value of these parameters. L nfortunatelv, the 
study of both laminar and turbulent flows have far outpaced the understanding of the 
transition mechanisms which link them. The transition phenomenon is believed to 
consist of three stages . 1 The first of these is the transfer of external disturbances into 
the viscous layer, which has become known as the receptivity problem. Next is the 
linear growth of these disturbances, and the final stage is the nonlinear breakdown 
to turbulence. The linear stability theory is well established and was first verified 
experimentally by Schubauer and Skramstad . 2 Extensions of the original theory to 
compressible and three dimensional flows have been made and efficient numerical 
schemes have been developed . 3 Advances in describing the third, nonlinear stage 
have more recently been made due largely to advances in computational capabilities. 
Direct numerical simulations of the transition process have been achieved by several 
authors.' 1, 6 

The majority of past numerical and theoretical work in transition, however, has 
relied upon simplifying assumptions such that temporal and spatial disturbances are 
considered separately. Furthermore, spectral methods have been extensi\elv em- 
ployed. These spectral methods are known to have difficulties with discontinuities, 
such as those found in compressible flows and flows over complex geometries. The 
present work extends the capabilities ot numerical solution of transitional flows by us- 
ing algorithms in wide use in computational fluid dynamics (CFD) which are capable 
of solving compressible flowfields. Furthermore, both spatial and temporal simula- 
tion is employed. This is an extension of the work by Fenno, Streett and Hassan' in 
which the transitional flow over a flat plate was numerically simulated. Receptivity 



in that work was of the forced type, behind the leading edge. The more realistic 
situation of pressure gradient flow and naturai leading edge receptivity is addressed 
in the present work. Lin. Reed and Saric 5 solved the case of incompressible flow 
over a flat plate with elliptic leading edge, which is characterized by a small region 
of favorable pressure gradient in the leading edge region, lollowed by an extended 
region of adverse pressure gradient. Their receptivity was of the natural type, with 
the external disturbance provided by an acoustic wave. The work was extended by 
Buter 0 to include receptivity to freestream vorticity disturbances. Murdock used an 
incompressible spectral method to study the leading edge receptivity of flat plates to 
an acoustic wave. He later extended his work to include parabolic cylinders, and thus 
favorable pressure gradient. 11 The previous four references worked with the incom- 
pressible problem. The current interest, however, is within the compressible regime. 
The objective of this work is therefore to numerically simulate leading edge recep- 
tivity of a parabolic cylinder using a compressible numerical scheme. The Mach 
number considered here is 0.3, for comparisons and validation to the available body 
of incompressible results. 


Numerical Scheme 


The governing equations are the two-dimensional, unsteady, compressible 
Navier-Stokes equations. A finite volume formulation using second order accurate cen- 
tral differencing is employed. The method follows the approach of Jameson, Schmidt 
and Turkel 12 and integration is by a four stage, modified Runge-Kutta method. Third 
difference explicit dissipation is used for numerical stability. The dissipation is scaled 
by the local eigenvalue of the equations. Following [12], the dissipation scaling pa- 
rameter is 

( 1 ) 


h i+±.j 

At 


where /?._i , is the cell height at (i + ^.j) and Af is the integration time step. In one 
dimension, the time step can be written as 

h : 


At = A t. = 


‘i+to 


( 2 ) 


ju| + a 

where u is the velocity in the x direction, a is the speed of sound and, therefore. 
|n| + a is the maximum eigenvalue in the x direction. Thus the numerical dissipation 
is scaled to the local eigenvalue. However, in two dimensions. 

1 1.1 
At ~ Af, ~ A t v 


( 3 ) 
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Simplification of equation (3) leads to 


(|u| + a) + (|u| + a) 

Therefore, the isotropic dissipation obtained by using equation (4) for the scaling 
term given in equation (1) leads to excessive dissipation. Jou et al. 13 have shown 
significant degradation of the laminar solution if excessive numerical dissipation, such 
as that described above, is involved. Improved models have been developed which 
use an anisotropic approach and are are only modestly more costly to solve than 
the standard scheme. 14, 15 Therefore an anisotropic dissipation model will also be 
investigated in the present study. This involves substitution of the directional A t x 
shown in equation |2) for A t in equation (1)- Likewise, the y direction dissipation 
is scaled by At y . This results in a slight increase in required computer resources 
because both direct ional time steps have to be either stored or recalculated at each 
grid location. It should be noted that the various forms of the time step mentioned 
above are for computation of the numerical dissipation only. The temporal integration 
time step uses the complete A t given by equation (4). 

The solid wall boundary uses no-slip and adiabatic boundary conditions. Ex- 
trapolation is applied at the outflow boundary, while the inflow boundary is handled 
using locally steady, one-dimensional characteristic theory. The positive half domain 
is solved and thus, symmetry boundary conditions are applied at the centerline. 

The first step in the solution procedure is to discretize the domain such that 
the relevant flow structures can be resolved. This is achieved through choice of an 
appropriate grid. The streamwise resolution must be fine enough to resolve the flow 
gradients in the mean flow near the leading edge. Streamwise grid spacing in this 
region is one tenth of a nose radius. This spacing requirement was found by numerical 
experimentation. The grid is stretched to a spacing no larger than one eighth of a dis- 
turbance wavelength at the outflow in order to resolve the desired spatial disturbance. 
In the body normal direction, the traditional rule is to place at least ten points in 
the viscous boundarv layer. In the present method, however, a more rigorous require- 
ment is the need to resolve the unsteady disturbance profiles. These profiles are of the 
acoustic and instability waves. The acoustic profile can be estimated with the Stokes 
wave, while the instability profile is estimated by linear theory. Both of these profiles 
have strong gradients and inflection points in the near wall region. A minimum of 4 
to 5 points is placed before this inflection point. A further requirement is the need 
to resolve the strong mean flow gradient at the leading edge. Grid spacing normal to 
the leading edge is one hundredth of a nose radius. In the temporal direction, two 
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constraints must be met. These are the need to resolve the temporal frequency and 
the need to meet numerical stability requirements. The more restrictive ot these is 
the numerical stability. 

The next step is to compute a mean flow. For this part ot the calculation, tem- 
poral accuracy is not needed and the integration is used only to achieve a steady 
state solution. Therefore, convergence acceleration in the form oi local time stepping 
and residual smoothing is used to hasten the computation. 16 It is extremely impor- 
tant that, for the high Reynolds numbers employed, a solution with smooth second 
derivatives be obtained. 1 ' 

Once the mean flow is obtained, the unsteady forcing disturbance is determined. 
The disturbance takes the form 

Aq = Aq e' {ax -“ t] (5) 


where a is the disturbance wavenumber, uj is the frequency, and 



(6) 


where e 0 is the specified disturbance amplitude. Kquation (6) is obtained from the 
one dimensional, inviscid conservation of mass and momentum equations and the 
definition of the speed of sound. 

The fourth and final step involves superimposing the disturbance on the free- 
stream at the inflow boundary of the computation domain. Thus a disturbed property, 
q, is written as 

q = Qoo + -^q U) 

where Qoo represents the freest ream value. This property becomes the inflow value 
used in the characteristic boundary conditions. The governing equations are then 
integrated in a time accurate manner, allowing the disturbance to continually prop- 
agate into the flow. Once the disturbance has penetrated the entire field and moved 
through the outflow boundary, sampling of the desired properties is performed. 


Results and Discussion 

Steady Results 

Results are presented for a Mach number of 0.3 such that comparisons can be 
made with the linear disturbance theory and other numerical results in the incom- 
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pressible regime. The problem is modeled after that ol Murdock 11 and the nose radius 
Reynolds number is set at R n = i x 10 J . A value of k = 2 x 10~ 3 for the coerhcient 
of third difference explicit damping was found to be sufficient to suppress odd-even 
point decoupling of the solution. The isotropic dissipation defined using equation (4) 
in equation (2) was initially used. The outflow boundary was set to R s 5 x 10 . 
A first order outflow extrapolation boundary condition based on Taylor series was 
applied along the grid lines. The farfield boundary was placed at x = —4 L upstream 
of the leading edge and y = 5665 above the body at the outflow. The length of the 
body within the computational domain is denoted by L, and S is the boundary layer 
thickness at the outflow boundary. The grid size was 416x352, which results in a 
resolution of between 374 cells per expected TS wavelength at the nose to 8 cells per 
wavelength at the outflow. The normal resolution uses a 1.65 percent stretching from 
a minimum spacing of y =<1/35. which places 29 points in the boundary layer. The 
results of the steady state solution are shown in figure 1. The figure shows the surface 


vorticitv defined as 


= -O 


C 2 

S, 


+ r 


( 8 ) 


where 0 is the nondimensional vorticitv and (f , rj) represent the grid coordinates. This 
parameter is shown vs. the log of if, and compares well with the results of Davis. 18 


Unsteady Results 

The unsteady disturbance was placed on the farfield boundary conditions and 
the integration was continued. This part of the computation does not use the conver- 
gence acceleration methods which were used for the steady state computation. The 
nondimensional disturbance frequency was set to F — 56 x 10 -6 . This is defined as 


where is the frequency, v is the kinematic viscosity and 0 is the freestream velocity. 
The amplitude of the disturbance was set to e 0 = 2x 10 -2 . This is of the same order 
of magnitude used by Buter and Reed 19 , who studied receptivity flow over blunt 
bodies and found the results to be linear up to forcing amplitudes of 4.2x10 2 . They 
attributed this to the favorable pressure gradient damping effect. 

The solution is compared to the results of Murdock 11 in figure 2. The temporal 
Fourier amplitude at the input frequency is plotted against the streamwise Reynolds 
number at a constant distance from the solid boundary given by a Reynolds number 
of Rr = 196. Qualitative similarity is seen between the present method and Murdock’s 
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spectral method results. The amplitudes show an initial strong spike near the nose, 
followed by a decay and a relatively constant value over most of the body. Murdock 
maintains that the magnitude of the envelope about the mean is associated with the 
Tollmien-Schlichting wave. 10, 11 

The instability wave profiles are shown in figure 1. normalized by the edge value 
of the total fluctuation and plotted against j] g . the Blasius similarity parameter. The 
instabilitv fluctuation is denoted by a double prime to distinguish it from the total 
fluctuation. With the exception of the first station at R, = 1 x 10 4 , the instability 
profiles are seen to take the characteristic dual lobe shape predicted by linear theory. 
The overall shape is seen to change from approximately equal magnitudes in both 
lobes to a dominant magnitude in the lower lobe. At about R , = I x 10 , the upper 
lobe has been reduced to insignificance. This single lobe shape is then seen to persist 
downstream. The instability fluctuation is seen to grow from about 10 percent of 
the total fluctuation at R s = 1 x 10 4 to almost 30 percent at R s = 8 x 10 4 . For 
comparison, the results of linear theory are shown in figure 4. Linear theory predicts 
growth rates, but not absolute magnitude. Thus, the profiles are normalized by their 
peak value. The different profiles were obtained by using the mean flow profiles at 
different streamwise stations obtained from the present method to obtain solutions 
using the linear theory code described in [ 1 ( ] - The linear theory produced similar 
profiles, as expected, for stations aft of R 3 = 4x 10 4 . but produced a differing profile 
at R s = 2xl0 4 . This is most likely due to the assumption of parallel flow used by the 
method, which is not valid in the leading edge region. The peak value is at a height 
of tj b = 1 for the theoretical profiles, whereas the present method profiles reach their 
peak at 77 3 v 0.6. The highly two-dimensional leading edge flow is also the most likely 
reason for the nonconforming profile shape at R s = 1 x 10 calculated by the present 
method and shown in figure 3. Furthermore, there is some region at the leading edge 
where the energy transfer between the long and short wavelength modes is occurring, 
and the idea of the acoustic and instability waves existing as separate waveforms can 
not be defended. Therefore no attempt should be made to analyze them separately. 
Goldsteim 0 argues that in the neighborhood of the leading edge, defined as 

, U°) 

the flow cannot be described as a double-layer structure of a Stokes wave and eigen- 
solutions of the Orr-Sommerfeld equation. In this region, the long wavelength free 
stream cisturbance is undergoing a wavelength reduction process by the non-parallel 
mean flow effects until the wavelength matches the Tollmien-Schlichting wavelength. 
Beyond this region, the disturbance wavelength is constant. In terms of the Reynolds 
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number, and usins the flow parameters of the present work, the limiting length in 
equation (10) can be expressed as 

R s < <9(l.8 x L0 4 ) . (11) 

which is in close agreement with the point where the uncharacteristic disturbance 
profile shape shown in figure 3 gives way to the proper shape. 

The maximum instability fluctuation of the momentum is shown in figure 5. In 
the figure, the maximum amplitude is determined over one period of the disturbance in 
time and also over the body normal direction. This yields a streamwise distribution. 
Three zones can be distinguished. The initial zone is \ery near the leading edge 
and is characterized by a sharp drop in amplitude. The next zone extends to about 
R s — 7 x 10 4 and indicates a strong increase in amplitude. The final zone extends 
to the outflow and shows a relatively constant amplitude with a long wavelength 
oscillation. These three zones correspond to the characteristic profile shapes given 
in figure 3. The first zone is associated with the profile at R s = 1 x 10 4 where the 
concept of separating the acoustic and instability waves is most likely not valid. The 
second zone is associated with the characteristically shaped instability profiles for 
1 x 10 4 < R, < R s = 7x 10 4 . The final zone, extending over the rest of the domain, is 
identified by the single lobe profile shown in the last profile of figure 3 for R s = i xlO 4 . 
Although the distribution shown in figure 5 is useful for identifying the zones described 
above, the quantitative amplitude information it offers must be treated as suspect due 
to the somewhat crude wave separation technique that was necessarily used to obtain 
it. The Stokes wave, however, is the best approximation of the acoustic signature 
that has been determined. 

The next step in the analysis is to determine the wavenumber of the computed 
instabilitv wave. Fisure 6 shows the instantaneous momentum fluctuation vs. the 
streamwise direction at a normal distance defined by 77 B = 0.6. This distance from the 
wall is where the largest fluctuation is seen in the instability profiles of figure 3. The 
long wavelength acoustic fluctuation is clearly seen in the figure, while a shorter length 
oscillation of much smaller amplitude is seen in the leading edge region. The smaller 
amplitude can be estimated as being about an order of magnitude smaller than the 
acoustic amplitude. This is in agreement with the magnitude of the instability profiles 
of figure 3. This shorter lengthscale oscillation is seen to be strongly damped beyond 
about R 3 = 5x 10 4 . 

The data series length represented by figure 6 is necessarily truncated with re- 
spect to the number of wavelengths. The acoustic wavelength is approximately equal 
to the bodv length, thus onlv one wavelength is available for analysis. Furthermore, 
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figure 6 indicates the existence of the instability wave tor only a tew wavelengths. 
These restrictions cause difficulties in the ability of traditional Fourier analysis tech- 
niques to resolve wavenumbers and amplitudes." 1 It was theretore necessary to use 
a technique better suited to the analysis of truncated series, the Maximum Entropy 
Method (MEM). 22 ' 23 ' 24 Corke" 1 maintains that this method is best used to detect 
the existence of various modes, rather than extract quantitative data. Figure 7 shows 
the results of the MEM when applied to the data series defined by figure 6. The 
power spectrum clearly displays three peaks. The first peak is at a wavenumber very 
close to the input acoustic wavenumber. The last peak is very close to the value pre- 
dicted by linear stability theory, illustrating the existence of the higher wavenumber 
Tollmien-Schlichting wave. The middle peak is associated with the numerical window 
function used in the analysis and is not indicative of any physical process. 

The existence of the TS wave is thus illustrated by identifying its wavenumber, 
at least qualitatively as seen in figure 7, and its profile shape, seen in figure 3. Its 
amplitude can be deduced from figure 6 as approximately an order of magnitude 
smaller than the acoustic wave and strongly damped in the streamwise direction. 
There are several possible sources of this damping. Since the TS wave propagation 
was successfully simulated using the present method for flat plate flows', complete 
with amplification and decay in the proper regions, it is reasonable to investigate the 
fundamental differences between flat plate flow and parabolic cylinder flow for the 
cause of the present TS wave damping. The first difference is the surface curvature. 
It is known that a convex surface has a stabilizing effect on the stability of a boundary 
layer, as compared to the flat plate. 25 ’ 26 However, Schlichting 26 also notes that this 
influence is very small if the ratio of the boundary layer thickness to the radius of 
curvature is small, as in the present case. Reshotko. 2 ° in his review of boundary layer 
stability refers to early work of Liepmann 27 who observed that transition on convex 
surfaces occurs at about the same Reynolds number as on flat plates. Therefore, while 
the curvature may be contributing to the TS wave damping, it is not believed to be a 
strong effect, especially since the curvature is small away from the leading edge. An- 
other factor which must be mentioned is the presence of a favorable pressure gradient. 
This has a much stronger effect on the stability of the boundary layer." 5- 26, 28, 29 The 
exploitation of a favorable pressure gradient as a transition control mechanism has 
even been studied. 30 Figure S compares the pressure gradient of the parabolic cylinder 
and the flat plate mean flows of [7], As expected, the pressure gradient is seen to 
be much more favorable for the parabolic cylinder until approximately R s = 1 x 10 5 , 
beyond which the results are essentially the same. This is well beyond the point 
where the short wavelength oscillation is observed in figure 6. This is strong support, 
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therefore, tor the pressure gradient as a leading cause of the instability v. a\e damping 
observed. 

Dissipation Model Study 

The results for the parabolic cylinder reported above used the two dimen- 
sional time step defined by equation |4) to scale the numerical dissipation. The 
anisotropic correction discussed following equation (4) was applied to the dissipa- 
tion at this point. The steady state flow was recalculated and the temporally accu- 
rate integration was conducted to 10 disturbance periods using an input disturbance 
amplitude of e o = 2xl0 -2 , as in the previous case. The acoustic amplitude at the 
body is e e = 8.8xl0 -3 , or 44 percent of the input amplitude. Previously, this was 
e„ = 1.3 x 10 -3 . or 6.5 percent. Figure 9 shows the TS profile for the two formulations 
at R s =bx 10 4 . Only a small change is seen in the instability response which can eas- 
ily be a result of the subtle changes in the mean flow. Thus, the improved damping 
formulation allows more of the input disturbance amplitude to reach the body, but 
the TS response remains essentially unchanged. 

Nose Radius Study 

Hammerton and Kerschen, 31 utilizing asymptotic methods, have shown the level 
of receptivity of a parabolic cylinder to freestream acoustic waves to be dramatically 
increased for lower Strouhal number. The Strouhal number is defined as 

S = FR n (12) 

and represents the ratio of nose radius to characteristic disturbance lengthscale. 
Therefore, to improve the TS response, the body nose radius was decreased from 
R n = 1000 to R n = 243. The Strouhal number is thus reduced from 5 = 5.6 x 10 -2 to 
5= 1.36x 10 -2 . 

The results of the unsteady computation show significant differences as compared 
to the Rn = 1000 case. Figure 10 shows the computed instability profiles. The 
stream wise distributions are compared in figure 11. As before, the first profile station 
at R s = 1 x 10 4 shows no characteristic instability shape and falls in the first zone 
of sharp amplitude decay in figure 11. The next two profiles fall within the second 
zone and show an initial growth of amplitude followed by a decay. Previously, only a 
growth was seen. Profiles between i? s = 7xl0 4 and i? s = 9xl0 4 develop a small lower 
lobe and a dominant upper lobe and lie within the third zone in figure 11. This third 
zone is also characterized by amplitude growth followed by decay. The final profile 
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at R s = i.l x 10 5 shows another shift of magnitude from the upper to lower lobe and 
again begins an amplitude growth which continues through the outflow boundary. 
This repeating growth and decay cycle was not seen previously for the R n = 1000 
case. That case does, however, show the growth portion of the first cycle before 
settling into a persistent single lobe shape. 

To explain the unusual behavior for this case, the mean flow pressure gradient 
distribution is displayed in figure 12 along with the distributions for the previous 
R n = 1000 case, and for the flat plate case. The pressure gradient, while still favorable, 
is much less so for the present case. This is a very convincing argument for the 
importance of the pressure gradient in the stability of the boundary layer. Whereas 
before no discernible TS behavior was seen beyond R s = l x 10 4 , here the TS profiles 
continue their characteristic shape to the outflow of the domain at R s = 1.2x 10 5 . 

Frequency Study 

The previous case involved a change in both the pressure gradient and the 
Strouhal number. To determine which of these changes is dominant, the final study 
is of a higher input frequency, F = 230x 10“ 6 and uses the R n =243 geometry. This 
frequency and geometry were chosen in order to maintain the same Strouhal number, 
S = 5.6x 10 -2 , which was used previously with the larger nose radius and lower fre- 
quencv. Furthermore, a higher frequency will produce a less stable boundary layer, 
as predicted by the upstream shifting of the branch points by linear theory. Further- 
more, Reed. Lin and Saric 8 have numerically shown a greater TS wave amplitude for 
higher frequency. The amplitude distribution is shown in figure 13 and shows very 
similar trends as for the lower frequency, larger nose radius results. This includes the 
initial amplitude drop, followed by an amplitude increase with the TS profile devel- 
opment. and then followed by a relatively unchanging amplitude for the remainder of 
the domain. The difference between the cases is in the range of R s where the various 
behaviors are seen. Previously, the TS profiles were in the range 2xl0 4 <R S < TxlO 4 . 
whereas in this case the range is 6 x 10 3 < R t < 1.6 x 10 4 . Thus, at a given R„ the 
higher frequency yields higher TS amplitudes. Note, however, that in terms of the 
nose radius, as seen in figure 14, both cases are very similar with the TS profile zone 
falling between 20 and 70 nose radii aft of the leading edge. The case which used 
S = 1.36xl0~ 2 . continues to follow its unique trends under the Rn normalization, how- 
ever. This serves to underline the importance of the Strouhal number, rather than 
either the nose radius or frequency alone, as the important similarity parameter for 
leading edge receptivity studies. 
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Concluding Remarks 

The objective of the present study was met in that the compressible Navier- 
Stokes equations, using a finite volume formulation, have been successfully used to 
compute leading ecise receptivity of boundary layer o\er a parabolic cylinder. Natural 
leading edge receptivity was simulated and Tollmien-Schhchting waves were observed 
to develop in response to an acoustic disturbance, applied through the farfield bound- 
ary conditions. The Tollmien-Schlichting waves which developed compared qualita- 
tively with the results of linear disturbance theory. The wavenumber of the TS waves 
was determined using the Maximum Entropy Method (MEM) and compared well 
with the wavenumber predicted by linear stability theory. 

The TS waves which were observed in the simulation were characterized by a 
strong amplitude attenuation in the spatial direction and the source of the attenuation 
was investigated. An anisotropic formulation of the artificial dissipation model was 
implemented. This allowed more of the input disturbance amplitude to reach the 
body, but the amplitudes were found to scale linearly to the original solution. Thus, 
the numerical dissipation was not found to have a significant effect on the strong TS 
wave damping. 

A smaller nose radius parabolic cylinder was employed and significant changes 
were seen in the unsteady results. To determine whether the changes were due to the 
resulting change in Dressure gradient or to the interaction of the acoustic wave with 
the smaller nose radius, a higher frequency case was computed using the smaller nose 
radius cylinder. The frequency was chosen such that the Strouhal number, or ratio 
of nose radius to disturbance wavelength, was the same as for the lower frequency, 
larger nose radius case. The solution was found to compare similarly to the original 
lower frequency solution in terms of the TS wave damping. Thus, it is concluded 
that the parameter which affects the development ot the TS waves is the Strouhal 
number, rather than the frequency or nose radius alone. The damping of the TS wave 
for the lower Strouhal number computation was, however, still significant. Therefore, 
the stabilizing effect of the favorable pressure gradient continued to have a dominant 
effect on the development of the TS wave. 

The present method, utilizing the compressible Navier-Stokes equations and a 
finite volume scheme, has been shown to be successful in simulating leading edge 
receptivity. The method has been validated in the nearly incompressible regime and 
for relatively simple geometries. With this accomplished, the extension of the method 
to higher Mach number flows and complex geometries is now possible. This capability 
has not previously been available. 
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Figure 5: Maximum Streamwise Momentum Fluctuation Minus Stokes Wave vs. Rey- 
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Figure 7: Fourier Power Spectrum of Spatial Streamwise Momentum Fluctuation 
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Figure 10: Maximum Streamwise Momentum Fluctuation Profiles Minus Stokes Wave 
for R n = 243, F = 56 x 10" 6 . 0.1 x 10 s < R, < 1.1 x 10 5 

20 

GWG5NAL PAGE IS 
OF POOR QUALITY 



0.50 — 


R n = 243 
Rn = 1000 


// / ! 
P^max!^ e 0.2 5 I 


' \ / 
t 

J 


x / \ / 
7 v 


/ 


0 .00 


0 . 0 


0 . 6 

R s x 10" 5 


1 . 2 


Figure 11: Maximum Streamwise Momentum Fluctuation Minus Stokes Wave vs. 
Reynolds Number for Different Nose Radii, F = 56 x 10 -6 
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Figure 13: Maximum Streamwise Momentum Fluctuation Minus Stokes Wave vs 
Reynolds Number for Various Frequencies and Nose Radii 
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Abstract 

A transition model for low-speed flows has previously been developed by Young, et al. 1 which 
incorporates first-mode (Tollmien-Schlichting) disturbance information from linear stability 
theory. The present work extends the model to high-speed flows by incorporating the effects 
of second-mode disturbances. The transition model is incorporated into a Reynolds- Averaged 
Xavier- Stokes solver with a one-equation turbulence model. Results using a variable turbu- 
lent Prandtl number approach demonstrate that the current model accurately reproduces 
available experimental data for first and second-mode dominated transitional flows. The 
performance of the present model shows significant improvement over previous transition 

modeling attempts. 
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Introduction 


The study of high speed transition is very important tor me efficient design ot hypersonic 
vehicles. Hypersonic vehicles, such as the National Aerospace Plane i X ASP 1. encounter flows 
which may be transitional over a significant portion of the vemcle. Within this transitional 
region, design parameters such as skin friction, heat transier. etc., are rapidly increasing 
or are maximum. Due to the current limitations ot supercomputers, the most practical 
method for describing such flows are methods based on the Reynolds- Averaged Xavier- 
Stokes (RANS) equations. Traditionally, transition models based on the Reynolds- Averaged 
equations have consisted primarily of the modification of existing turbulence models. The 
main objective of this work is the development of a transition model which incorporates 
information from linear stability theory and accurately reproduces available experimental 
data for transitional flows over a wide range oi Mach numbers. 

At Mach numbers less than about four, the transition process is dominated by first-mode 
disturbances. First-mode disturbances are vorticity disturbances which are characterized 
bv fluctuations in velocity. First-mode disturbances in incompressible flows are known as 
Tollmien-Schlichting waves. At Mach numbers greater than four, the transition process is 
dominated by second-mode disturbances. Second-mode disturbances are acoustical distur- 
bances which are characterized by large fluctuations in pressure and temperature, fluctua- 
tions which are much larger than the velocity fluctuations. 

Previous work by Young, et al. 1 has established a model for including the effects of first- 
mode disturbances into a transition model. The work by houng, et al. 1 dealt exclusively 
with incompressible flows. An objective of the present work is the extension of this model 
to account for Mach number effects. Additionally, the present work includes the effects of 
second-mode disturbances which dominate the transition process at hypersonic speeds. 


Experimental Data 

In order to determine the accuracy of models, it :s necessary to reproduce available experi- 
mental results. In the study of natural transition, flight data or data from quiet wind tunnels 
is necessary. The NASA-Langley Mach 3.5 Pilot Low-Disturbance tunnel has emerged as 
the facility capable of reproducing transitional data of the quality that matches flight data. 
Using this tunnel, Chen et alN presented recovery factor measurements on a o 3 half-angle 
cone and a flat plate. Three cases from these experiments are used for comparison with the 
present model. These are 

• Flow 3 - Case 5 Re = i.3e7/m x t = 0.0815m 

• Flow 3 - Case 6 Re = 5.39e7/m .r, = 0.1174m 

• Flow 3 - Case 7 Re = 3.35e7/m 


.r, = 0.2166m 


•:\IZ L ,■ = i-M.Tr,- = USA'. Mr..- = T 36 . </ = II. I 1 .: . The case designations are chosen 
:o correspond To the same cases presented in Ref. [3!. The reference quantities given are 
boundary- laver ease quantities. A.i oi the above data pertains to the cone experiments. 

m order to determine the eifect ot the seconti-moae disturbances model, experiments at 
Macn numbers greater than four are necessary. Kimmei et al. 1 have carried out Mach 8 
*ransition experiments over a (fared and non-ilared r half-angie cone. These experiments 
were carried out in the AEDC Tunnel B. The AEDC wind-tunnel is not a quiet tunnel 
in terms of freestream intensity levels. The ireestream turbulence intensity level in the 
experiments 4 was around 3.07c. however, only a small portion of this was in the second- 
mode frequence range (frequencies above 80-100 kHz). Based on this. Ref. [5] states, a 
conventional hypersonic wind tunnel can be quiet for second mode disturbances and noisy for 
first mode disturbances. The validity of this statement is still an open question. Regardless, 
at the present time, no second-mode dominated transition data obtained in quiet tunnels is 
available. Therefore, the present model is compared to these cases, which are given as 


• Mach 8 - Flared 

• Mach 8 - Xon-Flared 


M = 7.9S Re = 6562000/m 
.17 = 7.98 Re = 6562000/m 


Transition Model 

If r represents the fraction of the time the flow is turbulent, then the mean velocity, U m i s i 

r.„ = ru t + {i- U) 

In the above equation. Ut is the mean turbulent velocity and Ut is the mean non-turbulent 
or laminar velocity. Measurements by Kuan and Wang 3 showed that non-turbulent profiles 
are not Blasius profiles for flows over flat plates. Moreover, turbulent profiles are not the 
-raaitional fully developed profiles. At any given instant, the streamwise fluctuation in the 
Transitional region, u.. is given by. 

= u — Im 

If u is the non-turbulent velocity, tnen 

u r = l; — — U m = iif — Tilt — F/) 

and _ 


= Hi 


r 2 (U t -r,) 


-l2 


Similarly, if u is the turbulent velocity, then 


u'T = u‘ t 2 t { 1 — T) 2 (U t — Tf) 2 


( 2 ) 


(3) 
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i he expression given in equation :il is obtained 
l — l" ■ ot the time. Thereiore 


1' ul' the time ana that in equation t'2), 


- 


= r«i* - 1 1 — Du]' - rii - r )((:■- r f ) 


In general, one can write 


( U ;«') r = + (i - rhuith u + r(i - n a u, a u. 


(4) 


(5) 


The resulting expression is a result of three contributions: turbulent, non-turbulent. and 
a laree-eddv component. Traditionally, transition models have only included the turbulent 
contribution, terms resulting from non-turbulent fluctuations have not been included The 
terms A U, A U; are a result of large eddies. Their calculation requires specification of the 
turbulent and non-turbulent profiles in the transitional region. For this work, the large-eddy 

term has been neglected. , , , p., , 

The intermittencv factor. T. is obtained from the expression developed by Dhawan and 

Narasimha. b i.e. , ,v., 

T(j) = 1 — exp ( — 0.41 1^ 2 ) ( 6 ) 

with _ x-x t 

‘ “ A 

and A Is determined from the correlation 

Rex = 9 Re!? 


(7) 


x t is value of x at the beginning of transition. Because of this, this model requires the 
specification of a transition ’’point . 

Turbulent Contribution 

For a fullv turbulent flow, all that is needed to dose th e RAX S equation set is the turbulent 
Reynolds' stress. WR and the Reynolds heat flux. pu" } h" . By evaluating the turbulent 
Reynolds stress, the turbulent contribution to the general transitional stress is obtained. 

The turbulent contribution is the first term in Eq. 5. . 

The turbulent Reynolds stress is given in this work by Boussinesq s approximation which 

sives. . ... x .. 

(3) 


— PU’/U = jit 


dut 2 0u m 

doritr r ij dT 


m 


~ r^ijpk 


where p t is the eddy- viscosity and k is the turbulent kinetic energy per unit mass. 
In this work, a one-equation model is used which specifies the eddy-viscosity as 


Pt = A, p V k 


(9) 
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vnere - u.09 is a model constant ana i u is an aigebr?.:caily derived iengtii scale. Tnis 
iengtn scale is derived bv considering the near-wail behavior ot k and characteristics oi the 
ioe-iaw region in the turbulent veiocity boundary iaver. f’.e turbulent viscous length scale. 


> u . is given as 


( -,t = // 

1 — exp 

vM 
{ ' ( 

L 

^ j J 


(101 


where 

Ci = (0.3)"* k C 2 = 70.0 (11) 

and k is the Karman constant (0.411. The length scale C contains a damping factor which 
takes into account the presence of the laminar sublayer in the turbulent velocity profile. The 
form of the above length scale was originally derived by Gaffney. 


Reynolds Heat Flux 

In the energv equation, the quantity pu'jh" appears. This term is known as the Reynolds 
heat dux. This term is modeled by a gradient diffusion approximation, w mch gives 


pu'ih" 


„ - dT 

-GpCctP dx . 


( 12 ) 


where at is the turbulent diffusivity. The turbulent diffusivity can be expressed in terms of 
the eddy-viscosity and a turbulent Prandtl number. 


c*t 


Vj_ 

Pr t 


(13) 


The resulting expression for the Reynolds heat flux becomes 


pu'h" 


CpPt dT 

Pr> c Jx j 


(14) 


where C p is the specific heat at constant pressure, /q is the eddy viscosity, and Pr t is the 
turbulent Prandtl number. 

Traditionaiiv. the turbulent Pranatl number is taken as a constant. The value ot Pr t 
varies in the literature but generally is chosen between O.fi and 1.0. In this work, when the 
turbulent Prandtl number is chosen as a constant, the value is taken as 0.89. This was chosen 
to correspond to the value used in reierence [3]. 

There are physical consequences in choosing Pr t as a constant. By examining Eq. 14, it is 
clear that when Pr t is a constant, the damping factor used is the same factor which appears 
in ( u . This is equivalent to treating the turbulent thermal boundary layer with the same 
damping as the turbulent velocity boundary layer, even though they have different laminar 
sublayer characteristics. As results will illustrate, this choice has surprising consequences 
when considering quantities which are strong functions of the Prandtl number, such as the 
recovery factor. 
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Eddv-DifFusivity Model 

It is possible to model the Reynolds heat flux in such a way that the turbulent temperature 
boundary layer is damped differently than the turbulent velocity boundary layer. Based on 
the variable turbulent Prandtl number approacn of CebeciP the eudy-diffusivity moael is 
assumed to have the form _ 


n* = C 


A' t- 


1 r\ \ 


The temperature length scale, ir- is assumed similar in form to the viscous and dissipation 
length scales and is given as. 


( T — C3 y 



( 16 ) 


where 

C 3 = 3.12 
C 4 = 92.0 

The constants C 3 and C 4 were determined by comparing with the experimental turbulent 

Prandtl number measurements of Meir and Rotta. 9 

By defining a separate length scale for the temperature boundary layer, different damping 
is obtained. The present eddy-diffusivitv model implies a variable turbulent Prandtl number 
of the form. 

Pr, = ii (17) 

tr 

where is the viscous length scale. 


Non- Turbulent Fluctuations 

The second term in the general transitional stress. 

(\-T){uW). I 18 ) 

is a result of laminar or non-turbulent fluctuations. As stated previously, this non-turbulent 
contribution has traditionally not been considered in RANS-type approaches. In this work, 
the non-turbulent fluctuations considered are a result of first and second-mode disturbances. 

First-Mode Disturbances 

For moderate supersonic Mach numbers below approximately four, the dominant mode ot 
instability is the first mode. The dominant disturbance frequency at breakdown is well 
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■-redicred by the trequencv oi tfte nrst-mode disturDance bavins t;:e maximum ampiifica- 
• ;on rate. Using the work of Obremski et al.. 10 talker 11 nates that this frequency can be 
correlated by 

p = 3.2/fe.V (19) 

where U„ is the velocity at the ease ot ‘he boundary layer. Rex> is tne Reynolds number 
based on displacement thickness. ana is the irequency. 

The above expression is valid for low speed flow but can be extenaed. through the use 
of the reference temperature method, to high speed flow. It is known" that compressible 
formulas for skin friction and heat transfer have the same form as the incompressible formulas 
provided the physical properties oi the fluid are evaluated at some reference temperature 
between T w and T.. A possible definition of T * is 13 

T - = 1 - 0.032 M] + 0.58 ( R — l ) (20) 

-vhere the e subscript denotes boundary- layer edge quantities and w denotes wall values. By 
using this reference temperature. Walker s formula, i.e. Eq. 19, can be extended to higher 
Mach number flows by considering the definition of the unit Reynolds number as 

Re' = (21) 

with 

p = fx(T‘) (22) 

where T * is the reference temperature evaluated from Eq. 20. The density can be approxi- 
mated by notins 

R_ _ p 1 RT 1 
P m p~ R T~ 

and assuming that the pressure is approximately constant throughout the boundary layer. 


The dominant first-mode irequency correiation can then be rewritten as 


77 


= 3.2 ( Re’.) 


(24) 


with C m also evaluated using Re'. 

By assuming a scale of y l\ the flrst-mode disturbance frequency _• can be used to define 
a length scale. 


■ -.= = a 


a = 0.04 - 0.06 


(25) 
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It is assumed that the turbulent ana non-turbuient contributions to the stress tertsoi can be 
obtained by combining the two length scales as 

= ii - r-t„ - rc, 1 -6) 

where r is tne turbulent viscosity length scaie. Thus, when l = 0. i.e. before transition 


Second-Mode Disturbances 


Through the pioneering work of Mack. 14 - 15 the concept of second-mode disturbances was 
developed. As the boundary layer edge velocity increases above approximately M = 2.2, a 
region of the boundary layer becomes supersonic relative to the phase velocity. Mack showed 
that when such a situation exists, multiple solutions to the inviscid stability equations arise. 
These additional solutions are called the higher modes. The first of these modes is called 
the second-mode and is the most unstable of the higher modes, sometimes referred to as the 

Mack modes. 

At relatively low Mach numbers, boundary layer transition is dominated by lower fre- 
quency first mode disturbances. As the Mach number increases above four, the second-mode 
disturbances, which are higher frequency acoustical disturbances, dominate the transition 
process. From Mack 16 and Stetson and Kimmel. 5 the wave lengths of the most unstable 
second-mode disturbances have been shown to be about two boundary layer thicknesses. 

The frequency of the second mode disturbances can be written as 


where (A is the phase velocity and A is the wavelength. For hypersonic boundary layers, this 
leads to the relationship 

c Mi x .. )Q , 

\ X = — == l-S) 

V ft e T 

where c is a constant of proportionality. Combining with Eq. 2.. the frequency :s given as 


uJ = 


[ Ull (29) 

cM-x 

Again defining a velocity scaie as in Eq. 25. a second-mode contribution to the length scale 
can be defined as ^ 

f _ . 


W 


hich upon substitution of Eq. 29 above, yields 


■ ■ - [r. u - V 


kx 


(30) 


S 
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3y comparing with the Mach S stabiiitv experiment of Kimmei. 1 1 ' "he constant a was chosen 
as 0.23. Additionally, linear stabiiitv predicts the phase velocity to be aoout fl.'JU’.. Again, 
it is assumed that the second-mode contribution can be ootained by combining t>.u into the 
viscous length scale giving 

= .1 - n [(rs + Cm ) - r A (31) 

Results and Discussion 

Numerical Method 

The numerical method used to solve the equation set follows the method of Ref. [7] closely. A 
finite-volume method was used to evaluate the spatial terms in the Navier-Stokes equations. 
The inviscid fluxes were computed by an upwind approach based on Roe s flux difference 
splitting. MUSCL differencing was used along with a min-mod flux limiter. Second-order 
central differencing was used to evaluate the viscous fluxes. The solution was stepped in 
rime using a modified four-stage Runge-Kutta scheme until a steady-state was obtained. 
The determination of a steady-state was based on a density residual drop of five or more 
orders of magnitude. Due to the existence of a shock, grids of approximately 2 1 5x75 were 
needed to adequately resolve the shock and the very thin boundary- layer region. The steady- 
state was obtained wfith run-times on the average of 3-4 Cray Y-MP hours. 

Due to the very large Reynolds numbers, as high as 8x10', the thin-laver Navier-Stokes 
equations were solved. The thin-laver approximation neglects streamwise derivatives in the 
viscous terms due to their small magnitude relative to the normal derivatives. This is a valid 
approximation since at very large Reynolds numbers the viscous effects are restricted to a 
very small region near a solid boundary. Due to the cell aspect ratios necessary to resolve 
this small region, streamwise derivatives in the viscous terms have a negligible contribution. 

Mach 3.5 Cases 

The Mach 3.5 cases were used to determine the effectiveness of the transition model at pre- 
dicting transitional flows which are characterized by first-mode dominated transition pro- 
cessesf Previous work by Young et al. 1 has established the capability of the model for 
low-speed flat plate flows. The Mach 3.5 cases studied were the ones carried out on the 5" 
half-angle cone. A flat plate case was also studied but was not presented here. The cone 
experiments measured the recovery factor by determining the surface temperature under 
adiabatic conditions. 

The first set of results were obtained with the present model using a constant turbulent 
Pranatl number. The experiment data is from the Flow 3 - Case 5 experiment which is the 
first experiment of the Chen et al. 2 set. Fig. 1 compares the present model with the linear 
combination models of Xarasimha ana ONERA/CERT. Linear combination models assume 
the viscositv can be expressed as 

*- / ."V \ 
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The curve iabeied Xarasimha used the expression tor 1' -iven oy Lq. o wnnc Fie OX- 
ERA /CERT model is based on a transition function. «... wnich replaces T m 1 he above 
expression. This transition iunction. i,. is empirically derived ana is discussea m Kef. [3]. 
Fig. 1 is a plot of recoverv factor versus Reynolds number based on the distance along the 
surface. Fig. i demonstrates that the present model does a much better job of predicting 
the length of the transition zone as well as the peak in the recovery factor. The laminar 
region is also oredictea well. The Xarasimha and OXERA, CERT models were computed 
with a boundary laver coae and not a Xavier-Stokes formulation. It is believed that the 
boundary laver codes were started with a laminar profile just before the transition region 
and therefore the laminar region was never calculated. The recovery factor was assumed to 
be \/Pr and thus resulted in a straight line for the linear combination models. Both the 
linear combination models and the present constant Pr t model do a poor job of reproducing 


the decreasing trend of the recovery factor in the fully turbulent region. 

Fig. 2 compares the present results oi the constant turbulent Prandtl number model to 
the two-equation transition model of Wilcox. 1 ' The Wilcox model modifies the production 
and dissipation terms in the k-~ turbulence model to simulate transition. The Wilcox model 
does a better job of reproducing the peak in the recovery factor than the linear combination 
models. However, the present model with a constant Pr t does a slightly better job than the 
Wilcox model in calculating the recovery factor peak and does a much better job reproducing 
the transition extent. Additionally, the present method does a slightly better job in the fully 
turbulent region, although both fail to predict the downward trend in the recovery factor. 

The calculation of recovery factor is very sensitive to the choice of the turbulent Prandtl 
number. Fig. 3 is a plot of the current transition model with various choices of the constant 
Pr t . By varying the choice of Pr t by as little as 0.01. very large differences result in the 
calculation of the recovery factor. This is not too surprising since the recovery lactor is 
such a strong function of the Prandtl number. In Fig. 3, the solid line represents I he choice 
0 f p rt = 0.89 which is the choice used for the linear combination models as described in 


Ref. [3j. . .... 

Since the recovery factor is sensitive to the value of Pr-. a cioser look was taken at 

modeling the turbulent diifusivity. As stated previously, using a constant turbulent Prandtl 
number effectively treats the damping of the thermal boundary layer exacth the same as the 
damping of the velocity boundary layer. This has no physical basis, especially for high-speed 
flows where the thermal boundary layer becomes important due to the large temperature 
gradients in the viscous layer. Following the work of Cebeci. 3 the concept of a variable 
turbulent Prandtl number was used to define the current eddy-diffusivitv model ( Fq. 15). 
A typical profile of the variable turbulent Prandtl number is compared in Fig. & to the flat- 
plate measurements of Meir ana Rotta. 9 The results of this model for the Flow 3. Case 5 
experiment are shown in Fig. 4. As can be seen from the figure, the experimental data in 
the turbulent resion is aimost exactly reproduced. Even though the peak in the' recovery 
factor is slightly overpredicted, the transition extent is predicted very well. 

To determine if the constants used in the specification of the temperature length scale. 
Eq. 16. are general, the other experiments in the Mach 3.5 cases were computed. Fig. 5 is the 
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results of the present nioaei for T.e second case, r ;ovv 4 - * ase o. It is ciear mat the present 
moaei accurateiv reproduces the transition extent, the peak in the recover} factor, and the 
• irouient resion. The results snow rue extreme differences in the predictions ol the present 
metnod and the linear comDination models. A curious observation in the laminar region can 
■ ;je seen by examining Fig. 5. The experimental data is almost Hat in the laminar region while 
•he present method predicts a slight increase in the recovery factor. The experimental results 
mown here were obtained with the same test model as were the results in Fig. 4. However, 
different values were measured for the recovery iactor ior the same Reynolds number based 
on x. This seems to suggest either siight error in the measurements or the existence of 
some external phenomena. Fig. 6 is the results for the Flow 3 - Case i. Even though the 
transition extent is predicted well, poorer results are obtained in comparison with the other 
two cases. There is a possible explanation lor this. Each of the experiments were carried out 
with the same test model. A reasonaDle assumption is that the recovery factor measured 
should be the same lor each of the test cases lor the same Re T . given that laminar flow exists 
there. Fig. 7 demonstrates that each of the cases measured different values for the laminar 
recover v factor at the same Re-. As seen from the figure, the Flow 3 - Case 5 experiment 
measured around 0.837 for the recovery factor, while each of the succeeding cases measured 
values of about 0.840 and 0.S44 respectively. This is cleariv a discrepancy in the data. 
If the tests were not allowed to completely reach steady-state (adiabatic wall conditions), 
the measured value of the recovery factor would have been higher than the steady-state 
value. Since the computations assumed adiabatic conditions, and consistent values for the 
recoverv factor were computed in the laminar region, it is possible that the present method 
accurateiv predicts the steady-state value of the recovery factor for the laminar, transitional, 
and turbulent regions. 

The transitional results presented by Chen et al. 2 are well predicted by the present model. 
Moreover, the results of the current model illustrate the inadequacies of using a constant 
turbulent Prandtl number for cases that present recovery factor data. Evidently, the current 
•ransition model is the only model presently known, that correctly reproduces the trend in 
‘he turbulent region for these experiments." 11 

Mach 8 Cases 

The Mach S cases were usea to determine the effectiveness of the present model in predicting 
transitional flows which are characterized by second-mode dominated transition processes. 
The experiments were carried out by Kenneth Stetson and Roger Kimmel at the AEDC 
i Arnold Engineering Development Center l tunnel B and were reported in Ref. [4]. 

The first case considered is the flared cone case. The surface is flared to produce a con- 
stant adverse pressure gradient. Fig. !) presents the results for the present model without 
•he second-mode contribution giver, in Eq. 31. The results are presented as the dimensional 
surface heat flux versus distance along the surface. The laminar and the fully turbulent re- 
gion are slightly underpredicted for this case. Fig. 10 shows the same results for the present 
‘ransition model with the second-mode contribution. \\ it h this second-mode contribution. 
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-:;e iaminar cata is predictea correct;-/. However, the tuily-r urbiiierit region is siiii unaerpre- 
dctea. Since the surface was cooled, i: was determined that the modeled near-wail behavior 


>r the turbulent kinetic energy, given m Ref. 


Appendix Dj. needed modification. A factor 
, r ; ■ Jj— t.vas used in the present moaei to include effects ot wail temperature and to account 
for theVact that the near- wail behavior was originally modeled in Ref. [i] for adiabatic flows. 
The inclusion ot this factor results in the following expression tor the constant which appears 
in the turbulent viscous length scale i Eq. 10 h 


C 2 = wJf ( 33 ) 

Fig. 11 is the results of the present model with the near-wall modification. The experimental 
data is reproduced very well. The experimental data seems to imply that the heat flux 
in the turbulent region is either flat or is decreasing. This seems to be reinforced by the 
calculations in Ref. [17]. However, KimmeF points out that, “An adverse pressure gradient in 
incompressible flow causes boundary layer thickness to increase and heat transfer to decrease 
compared to zero-pressure gradient values, but the opposite trends occur in compressible 
flow. This is primarily because of changes in streamtube size and Mach number in pressure 
gradients in compressible flow. Consequently, wall shear and heat transfer decrease more 
slowly with x in an adverse pressure gradient than in zero pressure gradient. A strong 
enough adverse gradient causes heat transfer and wall shear to increase in the x-direction. 
Additionally. Kimmel 4 states that the error in the measurements could have been as high as 
—10%. It is therefore possible that the increase in the heat transfer shown by the present 
method in the turbulent region is accurate. 

The second case of the Mach 3 experiments was the non-flared 7° half-angle cone. Fig. 12 
presents the results of the current transitional model with the second-mode contribution as 
well as the near-wall modification. The results are presented as non-dimensional heat flux 
versus non-dimensional distance along the suriace. The transitional and turbulent regions 
are reproduced quite well. The discrepancy in tne laminar region is most likely experimental 
error since there is large scatter and no clear transition point. It is also possible that the 
laminar discrepancy is a resuit of the high levels of freestream turbulence. 

As the comparisons with the hypersonic cases demonstrate, the second-mode model per- 
forms quite well. 


Conclusions 

The present transition model is based on the fact that the stress in the transitional region can 
l/e expressed as a function of the turbulent and non-turbulent stresses. The non-turbulent 
stress is due to the presence of laminar fluctuations which are a result ot first and second- 
mode oscillations. This work has successfully formulated a method in which these laminar 
disturbances can be accounted for in the computation of a high-speed transitional flow. 
The previous work of Young et al. 1 has demonstrated the effectiveness of the first-mode 
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■ ontriburion for low-speed Hows, i"..e present work demonstrates that tne moael ran oe 
r-xtencieci to high-speed nows through the use oi the reierence temperature method. 

The present work successfully demonstrates that second-mode disturbances can also oe 
included into the transition model using a form similar to the nrst-mode disturbance model. 
The resulting transition model accurately reproduces available experimental data tor transi- 
tional flows. The moael has been shown to predict flrst-moae dominated transitional flows 
very accurately and performed better than previous modeling attempts. Additionally, the 
model has been shown to accurately predict hypersonic transitional flows which are charac- 
terized bv second-mode dominated transition. 
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Re x 


Figure 3. Flow 3 - Case a. Turbulent Prandtl Number Sensitivity 
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Figure 4. Flow 3 - Case 5. One-Equation Eddy Diffusivity Model Comparison 
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5. Flow 3 - Case 6, One-Equation Eddy Diffusivity Model Comparison 
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Figure 6. Flow 3 - Case 7. One-Equation Eddy Diffusivity Model Comparison 
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Figure 9. Mach 8 - Flared Cone. Present Method Without Second-Mode Model 
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Figure 10. Mach 8 - Flarea Cone. Present Method With Second- Mode Model 
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